* Multidimensional Identity Cleavages and Religious Discrimination *
* Replication DO File *

***********************************************************************************
* Table 1. Control Variables and Religious Discrimination, EGLS Panel Regressions *
***********************************************************************************

xtgls wsocdisx lxx v2x_polyarchy emajpct ln_zpop ln_gdp_pcap unempl_ilo emajchr emajmus i.year if !missing(erc)&!missing(rgc)&!missing(ric), corr(psar1) force
estimates store control1

xtgls wsocdisx lxx v2x_polyarchy emajpct ln_zpop ln_gdp_pcap unempl_ilo emajchr emajmus i.year if !missing(erxc)&!missing(rgxc)&!missing(rixc), corr(psar1) force
estimates store control2

xtgls mxx lxx v2x_polyarchy emajpct ln_zpop ln_gdp_pcap unempl_ilo emajchr emajmus i.year if !missing(erc)&!missing(rgc)&!missing(ric), corr(psar1) force
estimates store control3

xtgls mxx lxx v2x_polyarchy emajpct ln_zpop ln_gdp_pcap unempl_ilo emajchr emajmus i.year if !missing(erxc)&!missing(rgxc)&!missing(rixc), corr(psar1) force
estimates store control4

etable, estimates(control1 control2 control3 control4) keep(lxx v2x_polyarchy emajpct ln_zpop ln_gdp_pcap unempl_ilo emajchr emajmus) mstat(N) mstat(ll) note(Standard errors are in parentheses.) stars(0.10 "+" 0.05 "*" 0.01 "**" 0.001 "***") showstars showstarsnote center export(control.docx)

*********************************************************************************************
* Table 2. Cross-Cuttingness, Cross Fractionalization, and Religious Discrimination, EGLS Panel Regressions *
*********************************************************************************************

xtgls wsocdisx erc rgc ric lxx v2x_polyarchy emajpct ln_zpop ln_gdp_pcap unempl_ilo emajchr emajmus i.year, corr(psar1) force
estimates store all1

xtgls wsocdisx erxc rgxc rixc lxx v2x_polyarchy emajpct ln_zpop ln_gdp_pcap unempl_ilo emajchr emajmus i.year, corr(psar1) force
estimates store all2

xtgls mxx erc rgc ric lxx v2x_polyarchy emajpct ln_zpop ln_gdp_pcap unempl_ilo emajchr emajmus i.year, corr(psar1) force
estimates store all3

xtgls mxx erxc rgxc rixc lxx v2x_polyarchy emajpct ln_zpop ln_gdp_pcap unempl_ilo emajchr emajmus i.year, corr(psar1) force
estimates store all4

etable, estimates(all1 all2 all3 all4) keep(erc erxc rgc rgxc ric rixc lxx v2x_polyarchy emajpct ln_zpop ln_gdp_pcap unempl_ilo emajchr emajmus l.wsocdisx) mstat(N) mstat(ll) note(Standard errors are in parentheses.) stars(0.10 "+" 0.05 "*" 0.01 "**" 0.001 "***") showstars showstarsnote center export(all.docx)

*****************************
* Figure A. Point Estimates *
*****************************

coefplot (all1, keep(erc rgc ric)) (all2, keep(erxc rgxc rixc)), bylabel(SRD) || (all3, keep(erc rgc ric)) (all4, keep(erxc rgxc rixc)), bylabel(GRD) byopts(col(1)) xline(0) yline(3.5)

******************************************************
* Figure 1. Cross-Fractionalization Measures and SRD *
******************************************************

quietly xtgls wsocdisx erxc rgxc rixc lxx v2x_polyarchy emajpct ln_zpop ln_gdp_pcap unempl_ilo emajchr emajmus i.year, corr(psar1) force

margins, at(erxc=(0(.1)1)) post

estimates store erxc_srd

quietly xtgls wsocdisx erxc rgxc rixc lxx v2x_polyarchy emajpct ln_zpop ln_gdp_pcap unempl_ilo emajchr emajmus i.year, corr(psar1) force

margins, at(rgxc=(0(.1)1)) post

estimates store rgxc_srd

quietly xtgls wsocdisx erxc rgxc rixc lxx v2x_polyarchy emajpct ln_zpop ln_gdp_pcap unempl_ilo emajchr emajmus i.year, corr(psar1) force

margins, at(rixc=(0(.1)1)) post

estimates store rixc_srd

coefplot erxc_srd, bylabel(Religious-Ethnic CF) || rgxc_srd, bylabel(Religious-Geographic CF) || rixc_srd, bylabel(Religious-Economic CF) ytitle(SRD Index) byopts(row(1)) vertical recast(line) yline(0) lwidth(*2) ciopts(recast(rarea)) graphregion(color(white)) bgcolor(white) xlabel( 1 "0" 2 "0.1" 3 "0.2" 4 "0.3" 5 "0.4" 6 "0.5" 7 "0.6" 8 "0.7" 9 "0.8" 10 "0.9" 11 "1") legend(position(6) cols(3))

******************************************************
* Figure 2. Cross-Fractionalization Measures and GRD *
******************************************************

quietly xtgls mxx erxc rgxc rixc lxx v2x_polyarchy emajpct ln_zpop ln_gdp_pcap unempl_ilo emajchr emajmus i.year, corr(psar1) force

margins, at(erxc=(0(.1)1)) post

estimates store erxc_mxx

quietly xtgls mxx erxc rgxc rixc lxx v2x_polyarchy emajpct ln_zpop ln_gdp_pcap unempl_ilo emajchr emajmus i.year, corr(psar1) force

margins, at(rgxc=(0(.1)1)) post

estimates store rgxc_mxx

quietly xtgls mxx erxc rgxc rixc lxx v2x_polyarchy emajpct ln_zpop ln_gdp_pcap unempl_ilo emajchr emajmus i.year, corr(psar1) force

margins, at(rixc=(0(.1)1)) post

estimates store rixc_mxx

coefplot erxc_mxx, bylabel(Religious-Ethnic CF) || rgxc_mxx, bylabel(Religious-Geographic CF) || rixc_mxx, bylabel(Religious-Economic CF) ytitle(GRD Index) byopts(row(1)) vertical recast(line) yline(0) xline(1.5) lwidth(*2) ciopts(recast(rarea)) graphregion(color(white)) bgcolor(white) xlabel( 1 "0" 2 "0.1" 3 "0.2" 4 "0.3" 5 "0.4" 6 "0.5" 7 "0.6" 8 "0.7" 9 "0.8" 10 "0.9" 11 "1") legend(position(6) cols(3))

*********************
* Robustness Checks *
*********************

***********************************************
* Table A1. Summary of Descriptive Statistics *
***********************************************

sum wsocdisx mxx erc erxc rgc rgxc ric rixc lxx v2x_polyarchy emajpct ln_zpop ln_gdp_pcap unempl_ilo emajchr emajmus year if !missing(erc)&!missing(rgc)&!missing(ric)

*************************************************
* Figure A1. Cross-Cuttingness Measures and SRD *
*************************************************

quietly xtgls wsocdisx erc rgc ric lxx v2x_polyarchy emajpct ln_zpop ln_gdp_pcap unempl_ilo emajchr emajmus i.year, corr(psar1) force

margins, at(erc=(0(.1)1)) post

estimates store erc_srd

quietly xtgls wsocdisx erc rgc ric lxx v2x_polyarchy emajpct ln_zpop ln_gdp_pcap unempl_ilo emajchr emajmus i.year, corr(psar1) force

margins, at(rgc=(0(.1)1)) post

estimates store rgc_srd

quietly xtgls wsocdisx erc rgc ric lxx v2x_polyarchy emajpct ln_zpop ln_gdp_pcap unempl_ilo emajchr emajmus i.year, corr(psar1) force

margins, at(ric=(0(.1)1)) post

estimates store ric_srd

coefplot erc_srd, bylabel(Religious-Ethnic CC) || rgc_srd, bylabel(Religious-Geographic CC) || ric_srd, bylabel(Religious-Economic CC) ytitle(SRD Index) byopts(row(1)) vertical recast(line) xline(7) yline(0) lwidth(*2) ciopts(recast(rarea)) graphregion(color(white)) bgcolor(white) xlabel( 1 "0" 2 "0.1" 3 "0.2" 4 "0.3" 5 "0.4" 6 "0.5" 7 "0.6" 8 "0.7" 9 "0.8" 10 "0.9" 11 "1") legend(position(6) cols(3))

*************************************************
* Figure A2. Cross-Cuttingness Measures and GRD *
*************************************************

quietly xtgls mxx erc rgc ric lxx v2x_polyarchy emajpct ln_zpop ln_gdp_pcap unempl_ilo emajchr emajmus i.year, corr(psar1) force

margins, at(erc=(0(.1)1)) post

estimates store erc_mxx

quietly xtgls mxx erc rgc ric lxx v2x_polyarchy emajpct ln_zpop ln_gdp_pcap unempl_ilo emajchr emajmus i.year, corr(psar1) force

margins, at(rgc=(0(.1)1)) post

estimates store rgc_mxx

quietly xtgls mxx erc rgc ric lxx v2x_polyarchy emajpct ln_zpop ln_gdp_pcap unempl_ilo emajchr emajmus i.year, corr(psar1) force

margins, at(ric=(0(.1)1)) post

estimates store ric_mxx

coefplot erc_mxx, bylabel(Religious-Ethnic CC) || rgc_mxx, bylabel(Religious-Geographic CC) || ric_mxx, bylabel(Religious-Economic CC) ytitle(GRD Index) byopts(row(1)) vertical recast(line) yline(0) xline(3) lwidth(*2) ciopts(recast(rarea)) graphregion(color(white)) bgcolor(white) xlabel( 1 "0" 2 "0.1" 3 "0.2" 4 "0.3" 5 "0.4" 6 "0.5" 7 "0.6" 8 "0.7" 9 "0.8" 10 "0.9" 11 "1") legend(position(6) cols(3))

***********************************************
* Table A2. Control Variables and Religious Discrimination, PA Linear Panel Regressions with Stationary Correlations *
***********************************************

xtreg wsocdisx lxx v2x_polyarchy emajpct ln_zpop ln_gdp_pcap unempl_ilo emajchr emajmus i.year if !missing(erc)&!missing(rgc)&!missing(ric), pa corr(sta)
estimates store control1_xtreg_pa

xtreg wsocdisx lxx v2x_polyarchy emajpct ln_zpop ln_gdp_pcap unempl_ilo emajchr emajmus i.year if !missing(erxc)&!missing(rgxc)&!missing(rixc), pa corr(sta)
estimates store control2_xtreg_pa

xtreg mxx lxx v2x_polyarchy emajpct ln_zpop ln_gdp_pcap unempl_ilo emajchr emajmus i.year if !missing(erc)&!missing(rgc)&!missing(ric), pa corr(sta)
estimates store control3_xtreg_pa

xtreg mxx lxx v2x_polyarchy emajpct ln_zpop ln_gdp_pcap unempl_ilo emajchr emajmus i.year if !missing(erxc)&!missing(rgxc)&!missing(rixc), pa corr(sta)
estimates store control4_xtreg_pa

etable, estimates(control1_xtreg_pa control2_xtreg_pa control3_xtreg_pa control4_xtreg_pa) keep(lxx v2x_polyarchy emajpct ln_zpop ln_gdp_pcap unempl_ilo emajchr emajmus) mstat(N) mstat(ll) note(Standard errors are in parentheses.) stars(0.10 "+" 0.05 "*" 0.01 "**" 0.001 "***") showstars showstarsnote center export(control_xtreg_pa.docx)

***********************************************
* Table A3. Cross-Cuttingness, Cross Fractionalization, and Religious Discrimination, PA Linear Panel Regressions with Stationary Correlations *
***********************************************

xtreg wsocdisx erc rgc ric lxx v2x_polyarchy emajpct ln_zpop ln_gdp_pcap unempl_ilo emajchr emajmus i.year, pa corr(sta)
estimates store xtreg_pa_1

xtreg wsocdisx erxc rgxc rixc lxx v2x_polyarchy emajpct ln_zpop ln_gdp_pcap unempl_ilo emajchr emajmus i.year, pa corr(sta)
estimates store xtreg_pa_2

xtreg mxx erc rgc ric lxx v2x_polyarchy emajpct ln_zpop ln_gdp_pcap unempl_ilo emajchr emajmus i.year, pa corr(sta)
estimates store xtreg_pa_3

xtreg mxx erxc rgxc rixc lxx v2x_polyarchy emajpct ln_zpop ln_gdp_pcap unempl_ilo emajchr emajmus i.year, pa corr(sta)
estimates store xtreg_pa_4

etable, estimates(xtreg_pa_1 xtreg_pa_2 xtreg_pa_3 xtreg_pa_4) keep(erc erxc rgc rgxc ric rixc lxx v2x_polyarchy emajpct ln_zpop ln_gdp_pcap unempl_ilo emajchr emajmus l.wsocdisx i.year) mstat(N) mstat(ll) note(Standard errors are in parentheses.) stars(0.10 "+" 0.05 "*" 0.01 "**" 0.001 "***") showstars showstarsnote center export(xtreg_pa.docx)

*******************************************************
* Figure A3. Cross-Fractionalization Measures and SRD *
*******************************************************

quietly xtreg wsocdisx erxc rgxc rixc lxx v2x_polyarchy emajpct ln_zpop ln_gdp_pcap unempl_ilo emajchr emajmus i.year, pa corr(sta)

margins, at(erxc=(0(.1)1)) post

estimates store erxc_srd_xtreg_pa

quietly xtreg wsocdisx erxc rgxc rixc lxx v2x_polyarchy emajpct ln_zpop ln_gdp_pcap unempl_ilo emajchr emajmus i.year, pa corr(sta)

margins, at(rgxc=(0(.1)1)) post

estimates store rgxc_srd_xtreg_pa

quietly quietly xtreg wsocdisx erxc rgxc rixc lxx v2x_polyarchy emajpct ln_zpop ln_gdp_pcap unempl_ilo emajchr emajmus i.year, pa corr(sta)

margins, at(rixc=(0(.1)1)) post

estimates store rixc_srd_xtreg_pa

coefplot erxc_srd_xtreg_pa, bylabel(Religious-Ethnic CF) || rgxc_srd_xtreg_pa, bylabel(Religious-Geographic CF) || rixc_srd_xtreg_pa, bylabel(Religious-Economic CF) ytitle(SRD Index) byopts(row(1)) vertical recast(line) yline(0) xline(4.5) lwidth(*2) ciopts(recast(rarea)) graphregion(color(white)) bgcolor(white) xlabel( 1 "0" 2 "0.1" 3 "0.2" 4 "0.3" 5 "0.4" 6 "0.5" 7 "0.6" 8 "0.7" 9 "0.8" 10 "0.9" 11 "1") legend(position(6) cols(3))

*******************************************************
* Figure A4. Cross-Fractionalization Measures and GRD *
*******************************************************

quietly xtreg mxx erxc rgxc rixc lxx v2x_polyarchy emajpct ln_zpop ln_gdp_pcap unempl_ilo emajchr emajmus i.year, pa corr(sta)

margins, at(erxc=(0(.1)1)) post

estimates store erxc_grd_xtreg_pa

quietly xtreg mxx erxc rgxc rixc lxx v2x_polyarchy emajpct ln_zpop ln_gdp_pcap unempl_ilo emajchr emajmus i.year, pa corr(sta)

margins, at(rgxc=(0(.1)1)) post

estimates store rgxc_grd_xtreg_pa

quietly quietly xtreg mxx erxc rgxc rixc lxx v2x_polyarchy emajpct ln_zpop ln_gdp_pcap unempl_ilo emajchr emajmus i.year, pa corr(sta)

margins, at(rixc=(0(.1)1)) post

estimates store rixc_grd_xtreg_pa

coefplot erxc_grd_xtreg_pa, bylabel(Religious-Ethnic CF) || rgxc_grd_xtreg_pa, bylabel(Religious-Geographic CF) || rixc_grd_xtreg_pa, bylabel(Religious-Economic CF) ytitle(GRD Index) byopts(row(1)) vertical recast(line) yline(0) lwidth(*2) ciopts(recast(rarea)) graphregion(color(white)) bgcolor(white) xlabel( 1 "0" 2 "0.1" 3 "0.2" 4 "0.3" 5 "0.4" 6 "0.5" 7 "0.6" 8 "0.7" 9 "0.8" 10 "0.9" 11 "1") legend(position(6) cols(3))

*************************************************
* Figure A5. Cross-Cuttingness Measures and SRD *
*************************************************

quietly xtreg wsocdisx erc rgc ric lxx v2x_polyarchy emajpct ln_zpop ln_gdp_pcap unempl_ilo emajchr emajmus i.year, pa corr(sta)

margins, at(erc=(0(.1)1)) post

estimates store erc_srd_xtreg_pa

quietly xtreg wsocdisx erc rgc ric lxx v2x_polyarchy emajpct ln_zpop ln_gdp_pcap unempl_ilo emajchr emajmus i.year, pa corr(sta)

margins, at(rgc=(0(.1)1)) post

estimates store rgc_srd_xtreg_pa

quietly xtreg wsocdisx erc rgc ric lxx v2x_polyarchy emajpct ln_zpop ln_gdp_pcap unempl_ilo emajchr emajmus i.year, pa corr(sta)

margins, at(ric=(0(.1)1)) post

estimates store ric_srd_xtreg_pa

coefplot erc_srd_xtreg_pa, bylabel(Religious-Ethnic CC) || rgc_srd_xtreg_pa, bylabel(Religious-Geographic CC) || ric_srd_xtreg_pa, bylabel(Religious-Economic CC) ytitle(SRD Index) byopts(row(1)) vertical recast(line) xline(5.5) yline(0) lwidth(*2) ciopts(recast(rarea)) graphregion(color(white)) bgcolor(white) xlabel( 1 "0" 2 "0.1" 3 "0.2" 4 "0.3" 5 "0.4" 6 "0.5" 7 "0.6" 8 "0.7" 9 "0.8" 10 "0.9" 11 "1") legend(position(6) cols(3))

*************************************************
* Figure A6. Cross-Cuttingness Measures and GRD *
*************************************************

quietly xtreg mxx erc rgc ric lxx v2x_polyarchy emajpct ln_zpop ln_gdp_pcap unempl_ilo emajchr emajmus i.year, pa corr(sta)

margins, at(erc=(0(.1)1)) post

estimates store erc_grd_xtreg_pa

quietly xtreg mxx erc rgc ric lxx v2x_polyarchy emajpct ln_zpop ln_gdp_pcap unempl_ilo emajchr emajmus i.year, pa corr(sta)

margins, at(rgc=(0(.1)1)) post

estimates store rgc_grd_xtreg_pa

quietly xtreg mxx erc rgc ric lxx v2x_polyarchy emajpct ln_zpop ln_gdp_pcap unempl_ilo emajchr emajmus i.year, pa corr(sta)

margins, at(ric=(0(.1)1)) post

estimates store ric_grd_xtreg_pa

coefplot erc_grd_xtreg_pa, bylabel(Religious-Ethnic CC) || rgc_grd_xtreg_pa, bylabel(Religious-Geographic CC) || ric_grd_xtreg_pa, bylabel(Religious-Economic CC) ytitle(GRD Index) byopts(row(1)) vertical recast(line) yline(0) xline(3.5) lwidth(*2) ciopts(recast(rarea)) graphregion(color(white)) bgcolor(white) xlabel( 1 "0" 2 "0.1" 3 "0.2" 4 "0.3" 5 "0.4" 6 "0.5" 7 "0.6" 8 "0.7" 9 "0.8" 10 "0.9" 11 "1") legend(position(6) cols(3))
